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^O ■ Abstract 



We review the kinematics of a neutrino beam in the ideahzed case 
where the parent mesons momenta are parallel, but without any other 
approximation. This reveals several interesting features, in particular 

^~H, in the off-axis case, which are hidden by the approximations made in 

(— I ' a previous treatment. 
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1 Introduction. 

The kinematics of an off-axis neutrino beam in the approximation where 
the parent meson beam is divergenceless was posed as a problem to the 
student in these Archives a long time ago. However, the author of |T] seems 
not to have noticed that the problem is amenable to an exact treatment 
which reveals interesting aspects hidden by the approximations he uses. In 
this note, we show how to implement exactly the Lorentz transformation in 
the differential (angle, energy) distribution and we exhibit some features of 
the results. 

Among others, it is shown that a sizeable part of the flux at a given lab 
angle comes from neutrinos going backwards in the center of mass frame 
of the decaying meson. Given the approximations made in [1], it is not 
clear that they are taken into account in the results presented. How this 
might have a bearing on the calculations of neutrino fluxes being performed 
with more realistic divergent meson beams in the preparation of upcoming 
experiments is unknown to the present author. 

2 Conventions 

Starred variables refer to the c.o.m. frame (com) of the decaying meson. 
The neutrino energy in the lab is E and its angle with respect to the beam 
axis is 9 in the lab and 6* in the decaying meson frame. For illustrative 
purposes, we use a 30 GeV proton beam energy corresponding to the T2K 
experiment and a pion spectrum of the form {Ep — E^^)^ as done in [1] 

3 Neutrino generating decays 

vr^ ^ ^^ + n in) 

or 

K^TT + l {l) + Vl (Ul) 

with I = ij, or e and 

H ^ i^f, + e + i7e 

In the second and third cases, E* which is henceforth the neutrino com 
energy is not fixed by the conservation laws, but it makes little difference 



in what follows, except that the lab energy is not rigidly related to the lab 
angle and meson velocity. 

Also note that for jj, decays, the divergenceless hypothesis does not hold 
since /u's are themselves produced in meson decays. The forthcoming 
analysis should therefore be amended in this case. 

For the two-body case, we take n ^ n v as our paradygm: 
E^j_ + Ey = m-Tr and Pf^+ Pu = together with rrii, = imply that 

E^ = E* = "'^~"' and p* = E* = 29.8 MeV 



4 Lorentz transformation 

We take the decaying meson direction as polar axis in the following, and 
consider the meson beam as divergenceless. 

Consequently we shall not make use of the azimuthal degree of freedom and 
reduce the problem to two dimensions. 

The neutrino energy-momentum vector in com is therefore 
E*{1, cose*, sine*) 

Let 7 = ^ be the vr Lorentz factor and /3 = ^ its velocity for c = 1 

The lab neutrino energy-momentum vector is given by: 



E = jE* {1 + 13 cos 9*) (1) 

^ cos6' = J E* {cos e* + f3) (2) 

£; sin6l = E*sme* (3) 

We set p = E since our neutrinos are evidently considered as massless 
here. This is not an approximation when only these kinematics are consid- 
ered. To see, it is enough to remember that the neutrino Lorentz 7 in the vr 
rest frame is 3 10^ for a m = 1 eV/c^ neutrino, whilst the 7 of a 139 GeV 
vr is 10^ 

From this follows: 



cos 6* = ^ (4) 

smt^ = —, -— (5) 

7(l + /3cos0*) ^ ' 

All these equations are inversed by the substitution (3 -f-)- — /3 and 9 ■v^ 9* 

The relation for tan 9 following from dHandEI) has no interest. However, a 
complete picture of the relation 6 ^ 9* \s obtained by using tan a;/2 = iT^^ ^ 
which yields: 



tan 0/2= J I |tanr/2 

which shows that one angle is a monotonous function of the other on the 
complete [0,7r] interval, notwithstanding the well known headlight effect. 

The following useful relations should also be noted: 

(l-/3cos0)(l + /3cos0*)72 = 1 (6) 

E* 

which entails: E = — — - (7) 

7(1 — pcos9) 

cos 9-/3 

and the inverse of (HI): cos 9* = (8) 

1 — p cos 9 

5 Kinematical limits 

From equation 1^ one sees that sin 9 = ^- sin 9* and therefore 

E* 
sm 9 < 

- E 

which shows that a given neutrino energy E can only be found up to a 
maximum lab angle 

. / r.N F* 30 MeV 

t'max(-c;) = arcsm-g-w (9) 

which is small for most neutrinos. 



Conversely, at a given angle 6 from the (supposedly divergenceless) beam, 
the maximum neutrino energy is: 

EmUG) = —n (10) 

This bound is valid, of course, for not too small angles. Using ([7]), we 
can retrieve it more accurately by calculating ^ at fixed cos Q to get: 

dE E cose-j3 



O7 /37 1 — /3 cos 9 

It follows that the maximum lab energy at fixed cos 6 is obtained for 
/? = cos 6. This correspond to 7 = -^^ and we find, using again ([7]), that 
-Emax = -E*7 = ^j^ in conformity with (jlOp . These relations are illustrated 
on Fig. [1] 

So, although the bound seems to be independant of the beam, it makes 
sense only if there exist pions with E-,^ > mT^/svo.9 and (jlOp must be 
complemented by the restriction sin 9 > — - — 

The stationnarity of Ey['^Tr) will result in many neutrinos emitted by 
pions around that value of 7 to pile up at the same lab energy -Emax- Since 
they come from pions with 7's both below and above l/sin0, they are 
emitted both forwards and backwards in the pion frame. We will find the 
distributions of both components presently. 

For larger pion energies, neutrinos at angle 9 with respect to the axis 
are emitted more and more backwards in com with the result that their lab 
energy slightly decreases with E-,^ (see Fig. [T|) 

In an analogous manner, we can retrieve the bound ([9]) by using ([7]) to 
calculate the derivative §s^^ at fixed E. We find: 

d cos 9 1 cos 9 , 



and the minimum cos 9 at fixed lab energy is cos 9 = 13 so that again, the 
stationnarity wrt 7 will result in a piling up around cos 9 of neutrinos com- 
ing from pions with energies around m-,^/ sva.9. This is illustrated by the 
distributions (derived below) plotted on Fig. H] 
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E„„ at given iJ is obtained far 7, = 1 /sim5 (damned line — see text) 




Figure 1: Illustration of the -Emax(^) bound 



6 Probability distributions 

6.1 Fixed meson energy 

In this subsection, we consider that the decaying meson velocity is fixed. 
Therefore there is but one final state variable for a 2-body decay (not count- 
ing the azimuth) and the lab angular distribution is completely determined 
even in the case of 3-body decay of the parent meson. 



6.1.1 Angular distribution 

By differentiating ([5]) with respect to cos 9, we find 

acos^* _ 1-/32 

5cos0 ~ (l-/3cos6')2 



fill 



Since the meson is spinless it decays isotropically in com and the cos 9* 
distribution is flat on the interval [—1, 1] 

The lab distribution of the cosine for a given meson energy is thus: 



dP 1 1 - /32 

lip (12) 



acos^ 2{l-l3cos9)^ 

or, re-establishing the azimuthal degree of freedom 

^ = ' (13) 

dn 4^72(1-/3 cos 0)2 ^ ^ 

This distribution is valid without change for 3-body decays. 

6.1.2 Energy-angle distribution 

Equations ([T]) and ([4]) allow to calculate all the necessary partial derivatives: 



dE , ^ ^^. dcos9 

7(l + /3cosr) ^jT^ = 



dE* 'V ' ^ / QE 

dE ^^^ d cos 6 



dcos9* dcos9* 72(1 + /3 cos 0*)2 

from which we get the Jacobian: 

D{E*,cos9*) ,^ „ ,^, 1 

D[E^cos9) =7(l + /5-«n = ^(i_^,,,,) 

So assuming some energy distribution F{E*) for E* (independant of the 
angles by isotropy) we get: 

dP dP 1 1 

F{jE{l- 13 cos 9))- -(14) 



dEdn dE*dn*-f{l- l3cos9) " ' "^ 'M7r7(l - /3cos( 

For 2-body meson decays, there is only one variable {F is a Dirac 6) 
and this reduces to the angular distribution given in (jlSp after elimina- 
tion of 6; the result can alternatively be expressed as a distribution in energy. 



For 3-body meson decays, F must be found from the decay dynamics. 
From ([T|) and ([6]) -Emax for a given angle is given by (i_gcose') ^i^^ 

-^max — — ^ 2M obtained when the mass of the system recoihng 

against the neutrino in com is minimal. 

6.2 Introducing a meson spectrum 

We no longer consider that the energy of the parent meson is fixed. We 
assume an energy distribution, say g{E.,^) for the tt mesons in the lab. Note 
that we use 'vr' for definiteness. When it comes to 3-body decays, K would 
be more appropriate. 

Even assuming constant E* (2-body decay), the two final state variables 
E and cos 6 can now vary separately. To derive their distribution, we need 
the Jacobian of the transformation {Et^,cos6*) — t- {E,cos6) 



6.2.1 Transforming S^, cos 9* — > E, cos 9 

When performing this transformation, E* must be considered as a fixed 
parameter. In the case of 3-body decays, the {E, cos 9) distribution found 
below will be simply multiplied by the E* com distribution and integrated 
over this variable. 

It is much simpler to calculate the partial derivatives of the final variables 
with respect to E^^ and cos 9* than the contrary. All we need is again 
contained in the two relations ([T| and (JH) and the relation 

(15) 



d/3 1 






d'f (3'y^ 






The necessary derivatives are as follows: 






= —{1 + 1/ 13 cos 9*) 
dE^ m^ 


dE 
dcos9* 


= E*jf5 


9 cos 6* sin 6**2 1 


dcos9 
d cos 9* 


1 


dE^ {1 + /3 cos 9* fm^l3-f^ 


72(1-^/3 cos r)2 


Therefore we get: 







D{E ,cos9) _ E* COS0* 



D{E„, cos 9*) m^j'^ f3 1 + /3 cos 9* 
Using the relations established above, the inverse Jacobian reduces to: 

D{E^, cos 9*) m^p 1 



D{E ,cos9) E* {cos9-(3) 

which, not surprisingly, is singular when cos 9 = P, corresponding to 
cos 9* = 



6.2.2 Inverse transformation 

To get the final N(E, cos 9) distribution for given E* , we must invert the re- 
lation yielding (E, cos 9) as a function of (ii^vrj cos 9*). This is done as follows: 

Prom dll and ([5]) one easily deduce the following equations: 



(16) 
(17) 

r sin 9 
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1 — COS 9 cos 9* 






sin sin 0* 












(3 


= 


cos 9 - cos 9* 






1 — cos 9 cos 9* 




We further define: r 


— E 
= ~B 


- so 


that ([5]) reads sin 9 


Hence 


m 


and 


m 
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become: 
1- 








cos 9eVl — r"^ sin 
r sin 02 


L02 








/3 


= 


cos 

1- 


1 9 — e-y/l — T^ sin 


F 




cos 6'eVl — J"^ sin 


L02 



(18) 

(19) 

with e = ±1 according to whether cos 0* > or cos 0* < 0, that is, cos 9 > f3 
or cos 9 < /3 (cf. dH)); observe that ^T7\i implies that cos 9* is always smaller 
than cos 6*. Other constraints on /3 and 7 are automatically obeyed by the 
expressions ([16]) and ([Tj 



We also get from (fT9|) 

- ff sin ft f:^L/^ — r2 sin f)^ 

(20) 



cos 6 — P sin 0^e\/l — J'^ sin^^ 



P cos — e\/l — r^ sin 02 

This allows us to rewrite the Jacobian in terms of laboratory neutrino 
variables as: 

D{ETr, cos 9*) ttItt cos — eVl — r'^ sin 9'^ 
D{E , cos 9) ~^ sm9^Vl-r^sm9^ 

To see when each of the two a priori possible values of cos 9* apply, we 
observe first that ([1]) and ([7]) show that the physical range of r is 



Further, these two equations allow us to write cos 9 and cos 9* as func- 
tions of r and the boost parameters as follows: 

cos 9 



cos 9 



r7 — 1 



P^r 
The expression for cos 9* solves the e problem: 

• for r > 1 the two signs of cos 9* are possible since we can (if the beam 
permits !) have 7 < r and 7 > r The first case corresponds to a lower 
energy meson with the neutrino going forward in its rest frame, whilst 
the second means a more energetic meson needed to compensate the 
backward projection of the neutrino in this frame. 

• for r = 1 there is but one possibility left, corresponding to 

cos 9 = — cos 9* = 2^ 

PI 

• for r < 1, only the cos9* < solution exists. When r = I/7, cos 9 = 



J— 4 < T < y^l — /32 corresponds to neutrinos going backwards in 
the lab. 

This is where we diverge from the treatment of ^]. Equation (18) of 
this author seems to imply that only the cos 9* > contribution is taken 



10 



into account. However, cos6* = corresponds to 9 = arccos/3 which is 
generahy very smah, so that a sizeable number of almost forward neutrinos 
can be forgotten here. Fig. [2] strikingly illustrates this remark. 

Therefore, to get the joint distribution for (cos9,E) we must add the 
contributions of these two solutions when both apply. The joint distribution 
for E'tt and Q,* which is simply ^ \^^ is transformed to: 



ruy, 



dEdn ~ 47r E* sinO'^y/l - r'^sinO'^ 



g(?TT-7r7+)| cos^ — vl — r^ sin 0^1 + 5r(m^7_)| cos0 + v 1 — r^ sin 0-^1 

(22) 

where 7± means (fTSl) with e = ±1 and only the second term is to be kept if 
r < 1 . 

This must be multiplied by the probability density of E* and integrated 
over it to get the final distribution in {cos0,E). In the case of a 2-body 
decay, this results in the mere replacement E* — )• Eq in ([22]), where the 
value of Eq stems from the conservation laws. 

Figures [21 [3] and [H show the behaviour expected for constant E or con- 
stant cos 9. The last displays Jacobian peaks at 9 = aicsin E* / E correspond- 
ing to the limits explicited in ([9]) . The first shows the same phenomenon for 
fixed 9 with piling up at E = E*/sm9 in conformity with (jlOp The effect 
seems to disappear (Fig[3|) for 9 = .006. As explained after equation (fTO]l . 
the limit (in E) and the concomitant Jacobian peak exists only if pious of 
energy larger than 771,7^/ sin are present in the beam. With Ep = 30 GeV 
and the spectrum shape adopted here, this corresponds to 6* > 4.6 10^'^. 
However, the pion spectrum decreases sharply before this limit. 
Note that although the contribution of neutrinos going backwards in the 
pion rest frame dies off when going to smaller angles due to the scarcity of 
energetic pious, they concentrate at E*/sm9, reinforcing the peak which 
exists already in the distribution of neutrinos going forward in this frame. 

7 Conclusions 



Equation (|22p is the final result of this paper. It shows that there exist 
in general (for E > E*) two contributions to the spectrum of neutrinos 
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emitted at a given laboratory angle, provided the meson beam is energetic 
enough to bring the neutrino emitted backwards in the meson frame to 
the required forward angle. In fig. [2] this second contribution is displayed 
explicitely. 
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Figure 2: Distribution of E^, at fixed angle 
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Figure 3: Distribution of Ei, at fixed angle 
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Figure 4: Distribution of cos 9^ at fixed neutrino energy 



14 



